home *** CD-ROM | disk | FTP | other *** search
/ Fritz: All Fritz / All Fritz.zip / All Fritz / FILES / PROGMISC / PCSSP.LZH / PC-SSP.ZIP / STATMISC.ZIP / SRATE.FOR < prev    next >
Text File  |  1985-11-29  |  4KB  |  142 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C        SUBROUTINE SRATE
  5. C
  6. C        PURPOSE
  7. C           TO COMPUTE THE PROPORTION OF SUBJECTS SURVIVING, THE
  8. C           SURVIVAL RATES AND THE STANDARD ERRORS FOR SUCCESSIVELY
  9. C           REDUCED TIME PERIODS.  THE SURVIVAL RATE IS COMPUTED FOR
  10. C           EACH OF K PERIODS, WHERE K IS A CONSTANT TO BE SPECIFIED AND
  11. C           IS LESS THAN OR EQUAL TO N (WHERE N = TOTAL NUMBER OF
  12. C           PERIODS).
  13. C
  14. C        USAGE
  15. C           CALL SRATE (N,K,X,IE)
  16. C
  17. C        DESCRIPTION OF PARAMETERS
  18. C           N - THE TOTAL NUMBER OF PERIODS AFTER TREATMENT OR DIAGNOSIS
  19. C           K - THE SPECIFIED PERIOD UP TO WHICH SURVIVAL RATES ARE TO
  20. C               BE CALCULATED
  21. C           X - AN INPUT AND OUTPUT MATRIX (N X 9) CONTAINING THE
  22. C               FOLLOWING INFORMATION
  23. C               FOR INPUT--STORED IN THE N ROWS OF EACH COLUMN
  24. C                  COL 1 - NUMBER OF SUBJECTS ALIVE AT THE BEGINNING OF
  25. C                          PERIOD
  26. C                  COL 2 - NUMBER OF SUBJECTS WHICH DIED DURING THE
  27. C                          PERIOD
  28. C                  COL 3 - NUMBER OF SUBJECTS LOST TO FOLLOW-UP DURING
  29. C                          THE PERIOD
  30. C                  COL 4 - NUMBER OF SUBJECTS WITHDRAWN ALIVE DURING THE
  31. C                          PERIOD
  32. C                    FOR OUTPUT--STORED IN THE FIRST K ROWS OF EACH
  33. C                    COLUMN
  34. C                  COL 5 - EFFECTIVE NUMBER EXPOSED TO THE RISK OF DYING
  35. C                  COL 6 - PROPORTION WHO DIED DURING THE PERIOD
  36. C                  COL 7 - PROPORTION WHO SURVIVED DURING THE PERIOD
  37. C                  COL 8 - SURVIAL RATE
  38. C                  COL 9 - STANDARD ERROR OF THE SURVIVAL RATE
  39. C           IE- 1, IF K IS NOT IN THE CLOSED INTERVAL (0,N).
  40. C               2, IF THE NUMBER OF SUBJECTS ALIVE AT THE BEGINNING
  41. C               OF PERIOD I IS LESS THAN THE SUM OF THOSE WHICH DIED,
  42. C               WERE LOST, OR WERE WITHDRAWN DURING PERIOD I (I=1,...N)
  43. C               3, IF THE NUMBER OF SUBJECTS WHICH DIED, WERE LOST, OR
  44. C               WERE WITHDRAWN IN PERIOD I IS NOT EQUAL TO THE NUMBER
  45. C               ALIVE AT THE BEGINNING OF PERIOD I LESS THE NUMBER
  46. C               ALIVE AT THE BEGINNING OF PERIOD I + 1 (I=1,...N-1)
  47. C
  48. C        REMARKS
  49. C           IF THE SUBJECTS IN A GIVEN GROUP ARE ALL DIAGNOSED OR
  50. C           TREATED AT THE SAME TIME, THE CONSTANT K MAY BE SET EQUAL TO
  51. C           N .  IF THE SUBJECTS IN A GIVEN GOUP ENTER THE STUDY AT
  52. C           VARYING TIMES, K CAN BE NO GREATER THAN N-1.
  53. C
  54. C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  55. C           NONE
  56. C
  57. C        METHOD
  58. C           REFER TO S. J. CUTLER AND F. EDERER 'MAXIMUM UTILIZATION OF
  59. C           THE LIFE TABLE METHOD IN ANALYZING SURVIVAL', JOURNAL OF
  60. C           CHRONIC DISEASES, DECEMBER, 1958.  PP 699-712.
  61. C
  62. C     ..................................................................
  63. C
  64.       SUBROUTINE SRATE (N,K,X,IE)
  65. C
  66.       DIMENSION X(1)
  67. C
  68. C        INITIALIZATION AND ERROR CHECKING
  69. C
  70.       IE=0
  71.       NP4=4*N+1
  72.       NP9=NP4+NP4+N-2
  73.       DO 1 I=NP4,NP9
  74.     1 X(I)=0.0
  75.       IF (K) 2,2,3
  76.     2 IE=1
  77.       GO TO 45
  78.     3 IF(K-N) 4,4,2
  79.     4 DO 9 I=1,N
  80.       NP4=I+N
  81.       NP9=NP4+N
  82.       NP1=NP9+N
  83.       IF(INT(X(I)-X(NP4)-X(NP9)-X(NP1)+.01)) 5,6,6
  84.     5 IE=2
  85.       GO TO 45
  86.     6 IF(I-N) 7,9,9
  87.     7 IF (INT(X(I+1)-X(I)+X(NP4)+X(NP9)+X(NP1)+.01)) 8,9,8
  88.     8 IE=3
  89.         GO TO 45
  90.     9 CONTINUE
  91.    15 L1=0
  92.       L2=L1+N
  93.       L3=L2+N
  94.       L4=L3+N
  95.       L5=L4+N
  96.       L6=L5+N
  97.       L7=L6+N
  98.       L8=L7+N
  99.       L9=L8+N
  100.       LD=L2
  101.       LE=L5
  102.       LQ=L6
  103.       SUM=0.0
  104. C
  105.       DO 40 I=1,K
  106. C
  107. C        COMPUTE EFFECTIVE NUMBER EXPOSED TO RISK OF DYING
  108. C
  109.       L1=L1+1
  110.       L3=L3+1
  111.       L4=L4+1
  112.       L5=L5+1
  113.       X(L5)=X(L1)-(X(L3)+X(L4))/2.0
  114. C
  115. C        COMPUTE PROPORTION OF DYING
  116. C
  117.       L2=L2+1
  118.       L6=L6+1
  119.       X(L6)=X(L2)/X(L5)
  120. C
  121. C        COMPUTE PROPORTION OF SURVIVING
  122. C
  123.       L7=L7+1
  124.       X(L7)=1.0-X(L6)
  125. C
  126. C        COMPUTE SURVIVAL RATE
  127. C
  128.       L8=L8+1
  129.       IF (I-1) 20, 20, 25
  130.    20 X(L8)=X(L7)
  131.       GO TO 30
  132.    25 X(L8)=X(L8-1)*X(L7)
  133. C
  134. C        COMPUTE STANDARD ERROR OF SURVIVAL RATE
  135. C
  136.    30 L9=L9+1
  137.       SUM=SUM+X(L6)/(X(L5)-X(L2))
  138.    40 X(L9)=X(L8)*SQRT(SUM)
  139. C
  140.    45 RETURN
  141.       END
  142.